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Abstract. It is shown how some of the recent results of de Souza et aL can 
be generahzed to describe Hamiltonians whose eigenvalues are given as fe-generalized 
Fibonacci numbers. Here k is an arbitrary integer and the cases considered by de 
Souza et al. corespond to fc = 2. 



PACS numbers: 03.65.Fd, 02.10.De 

1. Introduction 

Curado and Rego-Monteiro introduced in j2] a new algebraic structure generalizing 
the Heisenberg algebra and containing also the g-deformed oscillator as a particular case. 
This algebra, called generalized Heisenberg algebra, depends on an analytical function / 
and the eigenvalues of the Hamiltonian are given by the one-step recurrence a^+i = 
/(a„). This structure has been used in different physical situations, see the references 
given in the recent paper [Ij. In the same paper [T] de Souza et al. introduced an 
extended two-step Heisenberg algebra having many interesting properties. In particular, 
they showed that in certain special cases the eigenvalues of the involved Hamiltonian 
are given by the well-known Fibonacci numbers, i.e., satisfy a two-step recurrence. It 
is the aim of the present note to show how one may introduce for arbitrary natural 
numbers k an extended k-step Heisenberg algebra which reduces for A; = 2 to the one 
discussed in pQ (and for A; = 1 to the one in [2J). In particular, the eigenvalues of the 
involved Hamiltonian are given in special cases by the /c-generalized Fibonacci numbers 
[3j . For the convenience of the reader we now recall the structure of the extended two- 
step Heisenberg algebra, using the notations of pQ. It is generated by the set of operators 
{H,a^ ,a, Js} where H = H^ is the Hamiltonian, a and - with a = (a^)^ - are the 
usual step operators and J3 = J3 is an additional operator. These operators satisfy the 
following relations: 

W = aHf{H) + J;), and aH = {f{H) + h)a, (1) 



ha^ = a^g{H), and aJ^ = g{H)a, (2) 
,at] = f{H)-H + J^, (3) 
H,,h] = 0. (4) 



The functions / and g are assumed to be analytical (in fact, the case discussed most 
thoroughly is the one where these functions are linear). In the Fock space representation 
one has the normalized vacuum state |0) defined by the relations 

a|0) =0, ii-jO) =ao|0), J3I0) = /3o|0), 
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where ao and /?o are real numbers. By some algebra one can show the following conse- 
quences 

a) \n) = Nn\n + I) , a\n) = Nn-i\n - I) , H\n) = an\n) , Jsjn) = /3n|n), 

(the first equation is the definition of the constants A'^) where the following relations 
hold true 

a„+i = f{an)+Pn, (5) 
Pn+l = 9{an), (6) 

^n+l = iVn + /(an+l)-an+l + /3n+l. (7) 

Combining @ and © shows the two-step dependence of the eigenvalues of the Hamil- 
tonian (i.e., the energy- levels): = /(«n) + gio^n-i)- In particular, choosing linear 

functions f{x) = rx and g{x) = sx (with r, s G M \ {0}) yields the recursion relation 
a„+i = ran + showing that the energy-levels comprise a generalized Fibonacci 

sequence (which reduces in the special case r = 1 = s to the usual Fibonacci sequence), 
thus providing a connection to Fibonacci-oscillators |4j. We may write the recursion 
relation of the linear C3jSG clS follows: 

'an+i\ _ f r I \ fa, 
Jn+i) ~\s 

Using j3n+i = san, we can write this also in terms of the {a„}„gN only (and where we 
switch the order of the vector components): 

T2{r,s) ^ =T2{r,s) ( (9) 



(8) 



.On+iJ \ s r J \ an J \ an J V ao 

where we have used the convention a_i = — . In p] the equation ([8]) is studied from a 
stability analysis and the fixed points are studied. For this the eigenvalues \±{r,s) = 
rifcVr^+4s q£ T2{r,s) are determined. However, it is also possible to use these to give 
an explicit expression for the Un due to ([9]). The result can also be derived by Binet's 
formula. Note that the case r = 1 = s of the classical Fibonacci numbers (where ao = 1 

and a_i = /3o = 0) gives rise to A-|-(1,1) = where A+(l,l) = is the golden 

ratio describing the limit lim„^oo "^''"^ of two consecutive elements of the sequence. 
One also has a relation between ([8]) and (generalized) Fibonacci chains via substitution 
rules [T]. In the case r = 1 = s one associates to ([8]) the substitution rule A — > AB 
and B ^ A, thus yielding the Fibonacci chain A AB ABA ABAAB ^ ■■■ 
whose length grows according to the Fibonacci numbers 1,2,3,5,.... In the case of 
arbitrary natural numbers r and s one has r + 1 different possible substitution rules, 
namely A ^ A - ■ ■ ABA ■ ■ ■ A (A appears r times and B appears at the i-th. place with 
1 < f < r + 1) and B ^ A---A {s times). 

In this note we show how some results of de Souza et al. [T] - corresponding to the 
case k = 2 - can be generalized to arbitrary natural numbers k > 3 where the integer 
k denotes the order of the recursion relation which the eigenvalues of the Hamiltonian 
satisfy. As recalled above, the case A; = 2 is associated to the well-known Fibonacci 
numbers which have appeared in numerous physical applications, many of them in the 
context of quasicrystals (through substitution rules - see, e.g., [3 [6] and the references 
therein as first steps into the vast literature) . The case of larger k is associated to the k- 
generalized Fibonacci numbers introduced by Miles [3] and discussed from many different 
points of view, see, e.g., [SEllHllinillllllllliaillllSllIlllIZlllSj. Up to now these 
numbers have appeared only in few physical applications (in particular in connection 
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with generalized exclusion statistics [191, '20] or the statistics of certain lattice animals 
which is closely connected to percolation problems |2H I22j). although they give also rise 
to interesting substition rules as will be discussed later on. For the convenience of the 
reader we have treated the case A; = 3 seperately from the general case, although this 
has led to some repetitions. 

2. The case A; = 3 

In this section we want to extend the structure discussed in the last section to an 
extended three-step Heisenberg algebra, i.e., to the case k = 3. We will accomplish this 
by introducing a further operator J4 with j| = J4 and a further analytical function 
h. The main idea is that the n-th eigenvalue of this new operator is connected to the 
(n — 2)-th eigenvalue of H (generalizing ^) and that the (n + l)-th eigenvalue of H - i.e., 
ctn+i - is given by f{an) and the sum of the n-th eigenvalues of J3 and J4 (generalizing 
(0)). Let us recall that ^ follows from ([2]) due to 

f3n\n) = Jsin) = -J—Jsa^n - 1) ® -J—a^ g{H)\n - 1) = g{an^i)\n). 

Thus, to obtain a connection between the n-th eigenvalue of J4 and the (n — 2)-th 
eigenvalue of H, we will have to consider what happens when (a^)^ is interchanged with 
J4. This motivates (jl2p . Similarly, to obtain the generalization of ([5]) we first recall that 
([5]) depends on ([1]) as follows 

an+i\n + 1) = -^W|n) ® -^a^fiH) + J3)\n) = {f{an) + /?„)[n). 

Thus, to realize the second part of the idea we have to generalize ([5]) by replacing on 
the right-hand side J3 by J3 + J4. This motivates (llOh . Of course, it is then necessary 
to close the algebra in a consistent fashion. After these motivations we now introduce 
the extended three-step Heisenberg algebra through the set of operators {H, ,a, J3, J4} 
satisfying the commutation relations (jl0|) -(|15|): 





= at (/(//) + J3 + J4), 


(10) 


jW 


= a^9{H), 


(11) 






(12) 


a, 




= f{H)-H + h + ,h, 


(13) 




] 


= 0,[i/,J4]=0, 


(14) 




] 


= 0. 


(15) 



(Here we have omitted the adjoint equations in (jl0p - (|12p .l In the Fock space represen- 
tation one has the normalized vacuum state |0) defined by the relations 

a|0) = 0, H\Q) = ao|0), J3IO) = /3o|0), J4IO) = 7o|0), 

where ao,/3o and 7o are real numbers. The first normalization constant is shown to be 
Nq = f{ao) — cto + /So + 7o (using (fT3l) ). By some algebra one can show the following 
consequences: 

at|n) = A''„jn + 1), a\n) = Nn^i\n — 1) , 
H\n) = an\n), J-sln) = f3n\n), J^ln) = jri\n). 
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It is now possible to show the following relations generalizing ([5]) and ([6]) 

Otn+l = /(On) + /3n + 7n, (16) 

(3n+i = g{an), (17) 
7„+i = h{a„-i), (18) 
Nn+i = A^^ + /(a„+i) -a„+i + /?„+! + 7„+i. (19) 
The new relation (jlSp is crucial for the following. The proof of it is analogous to (|17p : 

7n|n) = J4\n)= ^ Ma^f\n-2)^ ^ {a^fh{H)\n-2) 

jVn-liV„_2 iV„_iiV„_2 

(a))'^h{an-2)\n -2) = h{an--2)\n). 



Nn-l Nn-2 

Thus, 7„ = h[an-2)- Taking p6|) - (|18| ) together, we obtain the following relation 

an+i = /(an) + g{oLn-i) + h{an-2), (20) 

showing that the eigenvalues a„ of the Hamiltonian satisfy a three-step recurrence. 
Choosing linear functions f{x) = rx,g{x) = sx and h{x) = tx (with r, s,i G M \ {0}) 
yields the recursion relation 

ttn+i = ran + san-i + tan-2, (21) 

showing that the energy-levels comprise a generalized Tribonacci sequence (and in the 
special case r = s = t = 1 this reduces to the usual Tribonacci sequence). In analogy to 
the case k = 2 considered above, we can write the relations ()16p - (ll8p in the linear case 
as follows: 

' "n+i \ / r 1 1 \ / an \ 

Pn+l ]=[ S \ [ Pn ]. (22) 
7n+l / V i / V 7n / 

Note that we have used here 7n+i = ton-i together with /3„ = son-i to express 7n+i 
through As in the case k = 2 we can bring this into a form where only the {onlneN 
are involved (again, the order of the vector components is switched): 

1 \ / an-2 \ ( an~2 \ 

1 an-i =T3(r,s,t) a„„i . (23) 
t s r ) \ an ) \ an J 

Introducing for notational convenience the vectors a^~^^^ := (on-i a„,+i)* and using 
the conventions a^i = ^ and a_2 = we can iterate (I23p and obtain 

4"+^)=T3(r,.,i)"af . (24) 

It is the clear that the eigenvalues Xi{r, s,t), X2{r, s,t), X^i^r, s,t) of T3(r, s,t) will play 
an important role for the "dynamics" of the system (i.e., the solutions of (j2ip ) - as 
in the case k = 2 considered above; here we have assumed that the parameters r, s, t 
are chosen appropriately such that three eigenvalues exist. One also has a generalized 
Binet's formula (see, e.g., [31 (TTJ [TCI IE] ) which allows one to express the general solution 
On through the eigenvalues Xi{r,s,t) and the initial values aO)/3oi7o- However, note 
that a physical requirement will be that all energy-levels are greater than zero, i.e., 
En = q;„ > 0. This restricts the set of possible functions /, g, h, and, in the particular 
case of linear functions it restricts the coefficients r, s, t. For example, if r, s, t > than 
the energy- levels build a nondecreasing sequence, i.e., En+i > En- This seems to be a 
very natural condition (otherwise the counting seems to be very strange). 
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We can generalize the connection between ([8|) and (generalized) Fibonacci chains via 
substitution rules described at the end of the first section as follows. Starting from 
(j22p . we introduce the alphabet {A, B,C} and assume - in analogy to the case k = 2 - 
that r, s, t are natural numbers; in addition, we furthermore assume that 5 is a natural 
number, i.e., i is a multiple of s. Then there exist several substitution rules associated 
to ([22]): 

^ ^ A'Ml^'M2^'^ B^A', C^B-s, (25) 

where (^41,742) is a permutation of {B,C), < h < r for 1 < i < 3 and li + I2 + 
I3 = r (here we have denoted the concatenation of p letters A hj A^ = A - ■ ■ A {p 
times)). The number of substitution rules is given by (r + l)(r + 2) which is the 
number of different words of length r + 2 in the letters {A, S,C} where the letter 
A appears exactly r times. As an example, let r = 2,s = l,t = 2. A possible 
substitution rule is then given by A — s- ABAC,B — > yl, C ^ BB. This generates 
the following generalized Fibonacci chain: A ABAC ^ ABACAABACBB 
ABACAABACBBABACABACAABACBBAA ^ .... It is evident that the lengths 
grow very rapidly (1, 4, 11, 29, . . .). 



3. The generalization to arbitrary k > 2 



In this section we will generalize the situation of the preceeding sections to the case 
of arbitrary k > 2. To obtain a transparent formulation we will change the notation to a 
more convenient one. In addition to the Hamiltonian H and the step operators a and a"^ 
we have the operators Jj with jj = Ji for i = 2, . . . , k and associated analytic functions 
fi for 1 < i < k. It is clear from the explicit discussion of the case k = 3 that the main 
step consists in introducing new commuation relations for Jj and (a^)*~^. Thus, the set 
of operators of the extended k-step Heisenberg algebra consists of {H, a\a, J2, J3, ■ ■ ■ , Jk} 
and the commutation relations satisfied by these operators are given by ()26p -()30p: 



a, 



-A 



f,{H) + Y,jA 

\ 1=2 ) 



{a^y-^fi{H) for 2 < i < fc, 

k 



(26) 
(27) 
(28) 



[H, J.] 
[Jj, Jj] 



for 2 < i < /c, 

for 2 < I < j < /c. 



(29) 
(30) 



(Here the adjoint equations in ([26]) and (f27|) have been omitted.) Clearly, choosing 
A; = 3 reproduces the situation considered in the last section (if one also makes the 
replacement (J2, J3, /i, /a) ^ {Ja, Ji, f, 9, h)). In the Fock space representation one 
has the normalized vacuum state |0) defined by the relations 



a|0) = 0, H\0) 



a 



(0 



0) for 2 < i < k, 



(i) 

where Qq for 1 <i <k are real numbers. It follows that 



H\n) 



Nn\n + 1), a\n) = Nn^i\n - 1), 



al^^\n), Ji\n) 



a. 



n) for 2 < i < k. 
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The algebraic relations imply 



a. 



(1) 

n+l 

(i) 
n+l 



n+l 



i=2 



/.(«n-i+2) for 2 < i < /c, 



A^' + /i(«i'ii)-«i'ii+E«^li- 

i=2 



(31) 
(32) 
(33) 



As in the case k = 3 the relations (|32|) are crucial for the following and they are shown 
in the same way as (|18|) by application of (|27|) : 



n — i + l) = In ■ 

K{n, I) 



i + = fi{a, 



(1) 

n—i+l 



)\n). 



(Here we have denoted by K.{n,i) = Nn-iNn~2 • ■ ■ ^n-i+i the product of the constants 
appearing in the intermediate steps of the calculation.) Taking together (I3ip and (1320 
yields the relation 



«Si = /i(«i'^)+E/^(«n-m) 



(1) 



(34) 



1=2 



Let us now assume that the functions are linear, i.e., fi{x) = XiX for 1 < i < k (with 
Ai e R \ {0}). Then ^ reduces to 



a 



(1) 

n+l 



(1) 



(35) 



i.e., to the recursion relation of the k- generalized Fibonacci numbers [3l[TT]. In analogy 
to the cases A; = 2, 3 considered above we can write the relations (I3ip and ()32p in the 
linear case as matrix equation. For this we first observe that one has for 2 < i < k that 
ai+i = ^i<^n-i+2- Since on the right-hand side only expressions with index n should 



appear we have to reexpress a^_j^_2. 

„» _ A, (i-l) 



(i—1) 

However, one finds that ah. 



Ai-ia„_i+2' 



implying the sought-for relation al^!^_i = j-^On . To bring the recursion (|35p into 



the desired form we have to notice that a 



written as a„_j_i = AiOn + ah + 



(m) 
n 



Am a 



(1) 



so that this equation can be 



+ an ■ The resulting equation is 



(2T 



a, 



n+l 



a 



(3) 
n+l 



a, 



+1 
1) 



(4; 

n+l 



, (fc) 
V "n+l 



/ Ai 

A2 






1 


A2 






1 




A4 
A3 





1 


... 1 \ 


/ 


an 


\ 





... 




(2) 

ah 







... 




(3) 
ah 







... 




(4) 










ah' 









V 


(k) 
ah, 


J 



(36) 
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As in the cases k 
are involved: 



2 and A; = 3 we can bring this into a form where only the ^}neN 



a 



n-k+2 
(1) 

n— fc+3 



\ a. 



(1) 
n+1 











1 













As 



1 


A2 



\ 





1 

Ai ; 



a. 



(1) 

n-k+2 



V 



(37) 



Let us denote the matrix involved by T/c(A); here we abbreviate the coefficients Aj as a 
vector A = (Ai, . . . , Afc). 



(«i-fc+2 



Introducing for notational convenience the vectors 

(m + l) 
^fl 1 



•^n+i) using the convention a 

and obtain the generalization of 



(1) 



(n+1) 

we can iterate P7|) 



(38) 



It is the clear that the eigenvalues of Tfc(A) will play an important role for the "dynamics" 

of the system; however, the physical requirement of positive energies En = an ^ > win 
restrict the set of possible coefficients Aj as in the case k = 2> discussed above. One also 
has a generalized Binet's formula (see, e.g., [31 dH [HI [18]) which allows one to express 
the general solution through the eigenvalues of Tfc(A) and the initial values Let 
us mention that the matrix Tfc(A) is called in the mathematical literature (in particular 
in the case where all A, = 1) the k- generalized Fibonacci-matrix Qk [181 [TBI [T7] . 
If the coefficients Aj satisfy Aj > and Ai + • • • + A^ = 1 then Tfe(A) is a stochastic 
matrix and the recursion relation (I38p describes a Markov chain [14^ [T7] . If we denote 

I— I (k) 

the fc-generalized Fibonacci numbers introduced by Miles by FA , then one has in 



the case where all A 

ii) 



a, 



1 and where the initial values are given by Oq 
Q — for 2 < i < k the following relation 

(ai H h Cfc)! 



(1) 



E = n^' 



(1) 



n+k-1 



E 



0<a ]^ , . . ■ ,0-^ <n+fc — 1 
a I +2a2H \-kaj^=n 



ai! • • • Ofe! 



1 as well as 



(39) 



Here we have used the explicit expression for the A;-generalized Fibonacci numbers given 
by [3]: 

(ai H h ak)l 



E 

a-]^ +2a2 H — m — fc + 1 



oi! • • - afc! 



The case with arbitrary Aj as well as arbitrary initial values can be treated in a similar 
form using [llj . 

Clearly, it is again possible to associate substitution rules to (f36l) by introducing an 
alphabet {Ai,^2) • • • ,^k}- For this we have to assume that all Aj with 1 < i < k as 
well as all quotients Oj := with 1 < i < k — 1 are natural numbers. Note that 
this imphes A3 = q'3A2, A4 = q4^\3 = qiqi,X2 and in general Am = qmqm-i ■ ■ ■ q3^2- The 
appropriate substitution rules generalizing ([25]) are 



Ai > A-^ Ai-^A^ Ai^ ■ ■ ■ Ai^ 



.,A\\ A2 



^1 ' 



Ai 



Af_^ for 3 < i < /c 



(40) 



8 



MATTHIAS SCHORK 



where {ii,i2, ■ ■ ■ , ik-i) is a permutation of (2, 3, ... , k), < k < \i for 1 < i < k and 
h + h + • ■ ■ + l-k = ^1- The number of substitution rules is given by (Ai + l)(Ai + 
2) • • • (Ai + A; — 1) which is the number of different words of length Ai + fc — 1 in the 
letters {^i, . . . where the letter Ai appears exactly Ai times. 

4. Conclusions 

In this note we have shown how some results of de Souza et al. [T] can be generalized 
from the case k = 2 to arbitrary natural numbers A; > 3 where the integer k denotes 
the order of the recursion relation which the eigenvalues of the Hamiltonian satisfy. The 
case A; = 2 is associated to the well-known Fibonacci numbers which have appeared 
in numerous physical applications, many of them in the context of quasicrystals via 
substitution rules. The case of larger k is associated to the fc-generalized Fibonacci 
numbers introduced by Miles. It was pointed out that the case A; > 3 gives also rise to 
interesting substitution rules which might be interesting for the study of quasicrystals. 

References 

[1] de Souza J, Curado E M F and Rego-Monteiro M A 2006 J. Phys. A: Math. Gen. 39 10415 

[2] Curado E M F and Rego-Monteiro M A 2001 J. Phys. A: Math. Gen. 34 3253 

[3] Miles E P 1960 Amer. Math. Monthly 67 745 

[4] Arik M, Demircan E, Turgut T, Ekinci L and Mungan M 1991 Z. Phys. G: Particles and Fields 55 
89 

[5] Hida K 2004 Phys. Rev. Lett. 93 037205 

[6] Ilan R, Liberty E, Even-Dar Mandel S and Lifshitz R 2004 Ferroelectncs 305 15 

[7] Hoggatt V E and Bicknell M 1973 Fibonacci Q. 11 457 

[8] Fuller L E 1981 Fibonacci Q. 19 64 

[9] Philippou A N 1983 Fibonacci Q. 20 82 

[10] Philippou A N, Georghiou C and Philippou G N 1983 Internat. J. Math. Math. Sci. 6 545 

[11] Levesque C 1985 Fibonacci Q. 23 290 

[12] Dubeau F 1989 Fibonacci Q. 27 221 

[13] Lee G-Y and Lee S-G 1995 Fibonacci Q. 33 273 

[14] Mouline M and Rachidi M 1997 Rend. Semin. Mat. Messina, Ser. II 19 107 

[15] Dubeau F, Motta W, Rachidi M and Saeki O 1997 Fibonacci Q. 35 102 

[16] Lee G-Y, Lee S-G and Shin H-G 1997 Linear Algebra Appl. 251 73 

[17] Mouline M and Rachidi M 1999 Fibonacci Q. 37 34 

[18] Lee G-Y, Lee S-G, Kim J-S and Shin H-K 2001 Fibonacci Q. 39 158 

[19] Rachidi M, Saidi E H and Zerouaoui J 2003 Int. J. Mod. Phys. A 18 159 

[20] Schork M 2005 Int. J. Mod. Phys. A 20 4797 

[21] Turban L 2000 J. Phys. A: Math. Gen. 33 2587 

[22] Turban L 2001 |arXiv:cond-mat /0106595 1 



Alexanderstr. 76, 60489 Frankfurt, Germany 
E-mail address: mschork@member.ams.org 



